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Abstract 

We consider the uniform random d-regular graph on N vertices, with d € [Y“, for 

arbitrary a > 0. We prove that in the bulk of the spectrum the local eigenvalue correlation 
functions and the distribution of the gaps between consecutive eigenvalues coincide with 
those of the Gaussian Orthogonal Ensemble. 


1. Introduction and results 

1.1. Introduction. The universality of local eigenvalue statistics is one of the central questions 
in random matrix theory. Random matrix statistics are believed to apply to very general complex 
systems, including the zeros of the Riemann (^-function on the critical line. However, proofs of 
random matrix statistics have so far been limited mostly to matrix ensembles, with the notable 
exception [32]. There are two classes of matrix ensembles for which random matrix statistics 
have been established under very general conditions: invariant ensembles and ensembles with 
independent entries. For ensembles of random matrices that are invariant under the unitary or 
orthogonal group (invariant ensembles), much has been understood via the method of orthogonal 
polynomials (see e.g. [6ll 141138] 1 and, more recently, general results have been obtained by direct 
comparision of ensembles with different potentials muiQ]. For Wigner matrices and generalized 
Wigner matrices, whose entries are independent and typically nonzero, the universality problem 
has also essentially been solved completely [9]ll8p2f)lt22]l24p25ll29ll44| . For random sparse matrices 
with independent entries, significant progress has been made as well. In particular, for Erdos- 
Renyi graphs, in which each edge is chosen independently with probability p, random matrix 
statistics for both the bulk eigenvalues and the second largest eigenvalue of the adjacency matrix 
were established in m under the condition pN ^ jV^/s+o q- > Q. For the bulk 

eigenvalues, the lower bound on pN was recently extended to pN ^ A“ for any a > 0 in ISZj, and 
GOE statistics for the eigenvalue gaps was also established. Finally, aside from the approaches 
discussed above, supersymmetry has been used to obtain results on the local eigenvalue statistics 
for some special classes of distributions with independent entries (see e.g. m)- In addition, local 
random matrix statistics have been established by an analysis of transfer matrices (see e.g. m)- 
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In this paper we study random regular graphs, which are not invariant and do not have 
independent entries. We show that the eigenvalues of their adjacency matrices obey random 
matrix statistics in the bulk of the spectrum. The universality of local eigenvalue statistics 
for non-invariant matrix ensembles with correlated entries has recently been studied in a few 
other cases. In particular, after the appearance of this paper, GOE eigenvalue statistics were 
proved for the Laplacian matrix of sparse Erdos-Renyi graphs in [26], and for random matrices 
with certain short-range correlations in mm]; these results do not cover the hard constraints of 
random regular graphs. 


1.2. Main results. Let A be the adjacency matrix of the uniform random d-regular graph (RRG) 
on N vertices, i.e. a uniformly chosen symmetric matrix with entries in {0,1} such that all rows 
and columns have sum equal to d and all diagonal entries vanish. For d —>■ oo as —)• oo, it is 
known [3lll5tH5] that the eigenvalue density of (d — converges to the Wigner semicircle 

law whose density is g{x) ;= ^y^[4 — For d at least (logA^)^, three of the authors 

recently proved a local semicircle law for random regular graphs [3], giving precise estimates on 
the Green’s function and the eigenvalue density, down to spectral scales comparable with the 
typical eigenvalue spacing (up to a logarithmic correction). In this paper, we consider the local 
eigenvalue statistics of random regular graphs in the bulk of the spectrum. 

As the adjacency matrix of a d-regular graph, the matrix A has the trivial uniform eigenvector 
e := ...,!)* with eigenvalue d. We denote by Ai ^ ... ^ Atv-i the ordered nontrivial 

eigenvalues of (d — and by Errg the expectation with respect to the induced law on 

Xi ^ ^ Aat-i. By comparison, we denote by Eqoe the expectation with respect to the law 

of the ordered eigenvalues Ai ^ ... ^ Aat-i of the Gaussian Orthogonal Ensemble (GOE) on 
]^(Ar-i)x(iV-i)^ normalized so that the off-diagonal entries have variance N~^. 

The typical locations 7 * of the eigenvalues under the semicircle law are defined by 

^ = / Q{x)dx. (1.1) 

Theorem 1.1. Fix a > 0, and suppose that d G [N'^, Then, in the limit N —>■ 00 , the 

bulk gap statistics of the random d-regular graph coincide with those of the GOE. More precisely, 
for any fixed k > 0, n G N, and (f G C'“(]R"'), we have 

{Errg— ^goe) (l>{NQ{'yi){Xi — \i+i),... ,NQ{'yi){Xi — \i+n)) = o(l) (1.2) 

as N ^ 00 , uniformly in i G [kIV, (1 — «:)A^]. 


Next, let pj^ = p#,Ar denote the symmetrized joint law of the eigenvalues of the ensemble 
ff = RRG, GOE. The correlation functions are defined for n G |l,iV-ll by 

p^^\dXi,...,dXn) := p#(dAi,...,dA„,M^-i-). (1.3) 


Theorem 1.2. Fix a > 0 , and suppose that d G [N°^, Then, in the limit N —>■ 00, the 

locally averaged local correlation functions of the random d-regular graph coincide with those of 
the GOE. More precisely, fix a small enough constant c> 0, and define b = b]\f := Then 

for any fixed n G N, </> G C'^(M"'), and E G (—2,2) we have 
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For the GOE, the eigenvalue correlation functions are known explicitly; see e.g. |35j . Hence, 
the quantities for the GOE appearing on the left-hand sides of (II.2p and (jl.4l) can be computed 
explicitly. In fact, the eigenvalue gap distribution has only been computed in the sense of 
averages over the gap index; for the GUE, the computation for a hxed gap was performed 
in [i3] . 

The proofs of Theorems I1.1H1.2I follow the general three-step strategy developed in [20]- 
[22] : see e.g. [23] for a survey. In our setup, the strategy is formulated precisely in Section [2] 
The general idea is to study the convergence of eigenvalue statistics under Dyson Brownian 
motion (DBM) |16j . The three steps consist of (i) a local law providing precise estimates on 
the eigenvalue density down to the scale of individual eigenvalues, as well as the complete 
delocalization of the eigenvectors; (ii) the universality of the local eigenvalue statistics after the 
short time t = and (hi) effective approximation of the local eigenvalue statistics of the 

original matrix ensemble at t = 0 by the one evolved up to time t = 

In all previous instances of the three-step strategy outlined above, the independence of the 
matrix entries was crucial for steps (i) and (hi). Eor the random regular graph, a new approach is 
required for both of these steps, the last one of which is the main content of this paper. The local 
law for random regular graphs was recently established in |3], thus performing step (i). As for 
step (ii), the convergence of the local eigenvalue statistics under DBM with deterministic initial 
data was recently established in [33|, under the sole assumption that the eigenvalue density be 
bounded at the scale . Therefore the local semicircle law provides sufficient control on 

the eigenvalues so that using |3l[33| we can perform step (ii). 

Thus, the main difficulty is step (iii). There are several known methods for performing this 
step, including Lindeberg’s proof of the central limit theorem combined with higher moment 
matching conditions m, or the Green’s function comparison theorem [25|. Eor short times, a 
more direct method is to prove the stability of the eigenvalues under the DBM by analysing 
the dynamics of the individual matrix entries m- In all of these approaches, the independence 
of the matrix entries is used in an essential way. In contrast, the entries of random regular 
graphs are subject to hard constraints, and are therefore not independent. Tracking carefully 
the dependence of the matrix entries (using the methods from [3]), we find that the eigenvalue 
evolution is stable under a constrained DBM, for times t ^ Here, by stability, we mean 

that the changes in the local eigenvalue statistics are negligible. 

This stability can also be interpreted as follows: there is a class of reasonably well-behaved 
observables, which completely characterize the local bulk eigenvalue statistics, and whose time 
evolution under the constrained DBM can be well approximated by a switching dynamics of 
random regular graphs. We note that it has been proposed that, for random regular graphs, 
the dynamics provided by DBM should be replaced with a switching dynamics; see in partic¬ 
ular [HI]. However, obtaining rigorous results on the local eigenvalue statistics using only a 
switching dynamics is difficult, because the induced eigenvalue process is neither continuous nor 
autonomous [30]. Our strategy crucially relies on the fact that the eigenvalue process under 
DBM is continuous and satisfies an autonomous system of SDEs. 

Theorem ll.il and Theorem ll.2l hold also for sparse random matrices with independent entries; 
see m- We will use parts of that analysis which are applicable here. The main effort and novelty 
of this paper is in the control of eigenvalues under constrained DBM up to time t = using 

switchings. 

1.3. Further related results. Large regular graphs have been proposed as a testing ground for 
quantum chaos [42] . It is conjectured [7] that chaotic quantum systems (i.e. quantum systems 
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obtained by quantization of ergodic classical systems) exhibit random matrix statistics. Regular 
graphs are random matrices with a local structure, and as such a step in the direction of un¬ 
derstanding highly structured systems. It is believed that the eigenvalues of random d-regular 
graphs obey random matrix statistics for any d ^ 3. Indeed, there is numerical evidence that the 
local spectral statistics in the bulk of the spectrum are governed by those of the GOE [281137] . 
and further that the distribution of the appropriately rescaled second largest eigenvalue con¬ 
verges to the Tracy-Widom distribution of the GOE [36]. Our assumption d ^ for arbitrary 
a > 0 is purely technical, since some of the results used in our proof have only been established 
up to multiplicative errors of order N‘^ (with arbitrary c > 0). We believe that our results can be 
extended to d ^ (log with the same method. Furthermore, our results also extend to the 

other models of random regular graphs considered in [3|, such as the permutation model. Other 
results about the eigenvalue and eigenvector distribution of d-regular graphs on mesoscopic and 
macroscopic scales, with d —)• oo and with d fixed, are discussed in [3]. 

Our proof relies on switchings that leave the random regular graph invariant. Switchings 
of random regular graphs were introduced to obtain enumeration results in [M]; see m for a 
survey of subsequent developments. They are also used for simulation of random regular graphs; 
see e.g. m and references therein. Recently, switchings were used to bound the singularity 
probability of directed random regular graphs |12j . They also played an important role in our 
recent proof of the local semicircle law for random regular graphs [3|. 


Notation. We use a = 0{b) to mean that there exists an absolute constant C > 0 such that 
|a| ^ Cb, and a S> 6 to mean that a ^ Cb for some sufficiently large absolute constant C > 0. 
We use c for an arbitrarily small positive constant that may change from line to line. Moreover, 
we abbreviate |a, 6] := [a, b] n Z. We use the standard notations a Ab ■= min{a, b} and aV b ■= 
max{a, b}. Every quantity that is not explicitly a constant may depend on N, which we almost 
always omit from our notation. Throughout the paper, we tacitly assume N ^ 1. Unless 
otherwise stated, all sums of indices are over the set ll,Nl 


2. Strategy of proof 

Our goal is to prove that, in the bulk on the spectrum, the local eigenvalue statistics of A/y/d — 1 
are the same as those of the GOE. As mentioned in Section[Tl in order to show this, we interpolate 
between the RRG and the GOE using Dyson Brownian motion, or more precisely its Ornstein- 
Uhlenbeck version. 

2.1. Constrained Dyson Brownian motion. The adjacency matrix A of a regular graph is subject 
to the hard constraints that its rows and columns have fixed sum (i.e. it has the eigenvector 
e = ..., 1)*). Therefore, instead of the usual Dyson Brownian motion, we use Dyson 

Brownian motion constrained to the subspace of symmetric matrices whose row and column 
sums vanish. 

We begin with the notion of an Ornstein-Uhlenbeck process on a general finite-dimensional 
space. 

Definition 2.1. Let % be a real finite-dimensional Hilbert space. Let (f^jo he an orthonormal 
basis ofH. 
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(i) Let {wo,)a be i.i.d. standard normal random variables. Then we define the standard Gaus¬ 
sian measure on LL as W ■= '^^Wa^a- 

(ii) Let {ha)a be i.i.d. Ornstein-Uhlenbeck processes satisfying 

dha = dBa — -ha dt, 

where {Ba)a is a family of i.i.d. standard Brownian motions. Then we define the standard 
Ornstein-Uhlenbeck process on Ti as H{t) := ha{t) fa¬ 
it is easy to verify that the laws of W and the process H do not depend on the choice of the 
orthonormal basis (fa), and that the standard Gaussian measure is invariant under the standard 
Ornstein-Uhlenbeck process. We use these properties tacitly from now on. 

For example, let TL ■= {H £ : H = H*} be the Hilbert space of real symmetric N x N 

matrices with inner product 

{X,Y) := yTr(Xy). (2.1) 

Then the usual Wdimensional Dyson Brownian motion is the standard Ornstein-Uhlenbeck 
process H^t) on LL. More explicitly, H{t) is the Markov process satisfying the SDE 

dH = -^dB --Hdt, (2.2) 

y/N 2 ' ^ ^ 

where B{t) is Brownian motion on the space of x real symmetric matrices with quadratic 
covariation , Bki){t) = {dikSji + Sudjkfi. 

More intrinsically, given a finite-dimensional Hilbert space V, we denote the Hilbert space of 
symmetric linear maps on V with inner product ( 12 . 11 ) by TL{V). Then we define Dyson Brownian 
motion (DBM) on V to be the standard Ornstein-Uhlenbeck process on LLiV). With this point 
of view, the usual A^-dimensional DBM is the DBM on U = R^, and the constrained DBM is 
the DBM on U = e-*-. Note that the normalization N in (|2.1I) does not need to agree with the 
dimension of V, which is A^ — 1 for V = e-^. We make the convention to always normalize the 
inner product ()2.ip by N, no matter the dimension of V, and always denote the dimension of 
V by M. Finally, we denote the inner product on U by v • w for v, w G U. 

Definition 2.2 (Constrained DBM and GOE). The constrained DBM is the DBM on e-*-, 
i.e. the standard Ornstein-Uhlenbeck process onLLfe.^) with inner product (12.ip . The constrained 
GOE is the standard Gaussian measure on B(e-*-) with inner product (j2.ip . 

Thus, up to a change of basis, the constrained DBM is equivalent to the usual {N — 1)- 
dimensional DBM, with the minor difference of normalization by N rather than N — 1. However, 
since the definition of the d-regular graph is tied to the standard basis of R-^, it is frequently 
convenient to work with the constrained DBM in the standard basis of R'^. 

Next, in accordance with the decomposition R^ = e-*- 0 span(e), we have a canonical iso¬ 
morphism H H := B 0 0 from T-L{e^) to the set of matrices 

M := {BgR^""^ :B = B*,Be = 0}. (2.3) 

Throughout this paper, we tacitly identify H and H. 

We denote by C^{Ai) the space of functions F : Ad —>• C with continuous bounded derivatives 
up to order n. Sometimes it will be convenient to compute derivatives of functions F G 
in directions of R^^'^ that do not lie in Ad, which is made possible by the following convention. 
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Figure 2.1. A simple switching is given by replacing the solid edges by the dashed edges. 


Definition 2.3. Let P = I — ee* be the orthogonal projection from onto e-*-. We extend 
any function F G C”'(A4) to a C"'-function on through 

H ^ + H*)P^ , 

and denote this extended function also by F. Finally, for any F G C^{M.) and i,j G |1,A^], we 
use the abbreviation dijF{H) = -^—(Ff). 

From now on, we take W to be the constrained GOE and H = F[{t) to be the constrained 
DBM, with initial condition 


H{0) := / {A- dee*) G M. (2.4) 

Vd — 1 

Here A is the adjacency matrix of the random d-regular graph. In particular, the eigenvalues of 
H{0) as an element of F{e-^) are the rescaled nontrivial eigenvalues of A. 

2.2. Switchings. Simple switchings are an especially convenient generating set of Ad; they play 
a central role throughout this paper. For any i,j,k,l G |1,A^1 we define the simple switching 
efj' G A4 by 


^ij — ^ij ^kl ^ik ^jl where — dip6jq -|- diqdjp . (^-b) 

The action of a simple switching on an adjacency matrix, given by A i—)• A-\-^^j, amounts to 
adding the edges {i,j}, {k, 1} and removing the edges {i, k}, {j, /}; this is illustrated in Figure ITT] 
and made precise in (j3.14p below. In this section, the four vertices need not be distinct. 

Next, we define the abbreviations 

:= := = Tr(e^'5), (2.6) 

for all i,j, k, I G |1, A^]. Here dx denotes the directional derivative in the direction X. Explicitly, 
expressed in the standard basis on we have 

= 2{H,j + Hki - Hik - Hji ), (2.7) 

d^jFiH) = 2(5,, + dki - dik - dji)F{H ), (2.8) 

where F G C^(Ad). With these abbreviations, the generator of the constrained DBM can be 
expressed in terms of switchings as stated in the following proposition. 
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Proposition 2.4. The generator of the constrained DBM from Definition \2.S\ is 


L : = 


_ \ ^ /pikl\2 _ ^ \ ^ ffklai 

]\f3 32N‘^ ^ 


16A^3 




:ikl 

hj ' 




This means that for any F G C^(A4) we have 


d 


-E[F(/7(i))] = E[LF{H{t))] . 


(2.9) 


( 2 . 10 ) 


Proof. Let H(t) be the standard Ornstein-Uhlenbeck process from Definition 12.II on the space 
with inner product (| 2 .ip . As in the example (| 2 . 2 p , we obtain the quadratic covariation 

{Hij ,Hki){t) = j^{5ikdji + Sii6jk)t. (2-11) 

Next, let R G 0{N) satisfy Rejsf = e. Then, since the inner product (12.ip is invariant under 
orthogonal conjugations, we can express the constrained DBM as H{t) = R{H(t) 0 0)R*. We 
abbreviate H = H{t) and write for F G C^(A4), using ltd calculus, 

dEF{H) = -i^E[if,,(9i,F)(F)]dt + i E[(ai,-4ii")(if)d(i?i,-,Ffc,)] . 

i,j i,j,k,l 


By definition of R we have Rin = for all i, so that (j2.1ip yields 


N-l 


d{Hij,F[kl} — ^ ^ {fRimRjnRkmRln T RimRjnRknRlm) dt 

m,Ti=l 

= ^ik - (Sji - dt + ^(6u - (Sjk - Idt. 


N 


iv V A^y V N 

Thus, for any F G C^(A4) we have (I2.inh with 

^ ^ W <911 (<911 + dki - dii - djk) Hijd, 


N 




( 2 . 12 ) 


i,j,k,l 




Finally, using ^ij = ^ji = 0 for Lf G A4, we observe that L from (12.121) can be rewritten 
as (j2.9p . □ 

2.3. Outline of proof of Theorems II.1111.21 Theorems 11.1111.21 are an immediate consequence of 
the following two propositions. As in [3], we set 

iV 2 

D := dA^. (2.13) 

d-^ 

We always assume d G [A^“, which implies D ^ A^“. To state the two propositions 

concisely, we introduce the following definition. It will also be convenient in the proofs. 

Definition 2.5. Given H G M., we denote by Xi ^ ••• ^ Aw-i the eigenvalues of H\^a_. 
Consider two random matrix ensembles Hi and H 2 in Ai. Then we say that 
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(i) the bulk eigenvalue gap statistics of Hi and H 2 coincide if for any n £ N, 4> G 
and K > Q, we have 


[¥.Hi-^H 2 )<P{NQ{'yi){Xi-\i+i),...,NQ{'yi){\i-Xi+n)) = 0 ( 1 ) (2-14) 


as N ^ 00 , uniformly in i £ l^iV, (1 — «:)A^]; 

(a) the averaged bulk eigenvalue correlation functions of Hi and H 2 coincide if for any n £ N, 
(f> £ c > 0 small enough, and E £ (—2, 2), we have for b := 


^ [ dE' f f){xi,...,Xn)N^ {phI - PhI) {E' 

JE-b ilR" V 


dxi dXn \ 

^ Ng{E) ’ ■ ■ ■ ’ ^ Nq{E) ) 


= 0 ( 1 ), 

(2.15) 


where the correlation functions are defined as in ()1.3p . 

Moreover, we say that the bulk eigenvalue statistics of Hi and H 2 coincide if (i) and (ii) hold. 


Proposition 2.6. For any fixed <5 > 0 and t ^ N ^ , the bulk eigenvalue statistics of 

H{0) and H{t) coincide. 


Proposition 2.7. For any fixed <5 > 0 and t N the bulk eigenvalue statistics of Hft) 
and 77 ( 00 ) = W coincide. 


Propositions 12.6112.71 are proved in Section [5l As mentioned in Section (H our main effort 
and novelty is in proving Proposition 12.61 Proposition 12.71 is essentially a consequence of general 
results on universality of local eigenvalue statistics with small Gaussian component [33]. The 
local semicircle law of [3| is an important input in the proofs of both propositions. 

Proof of Theorem 11.21 The proof is immediate from Propositions [2Tll2.7[ with 5 ^ a/4. □ 


3. Switchings and short-time comparision 


The main result of this section is Proposition 13.II below. To state it, we introduce the following 
Sobolev-type seminorms, whereby the derivatives are taken in the directions of all switchings 


A := £R^^^ ■.i,j,k,l£ll,Nl}. (3.1) 

First, for r ^ 1, we define an L'’-seminorm on C^{Ai) through 

\\F\\r,t := (E|F(77(^))r)'/^ (3.2) 

Then, we extend this seminorm to include derivatives; for F £ C^(A4) we define 




r,t 


sup sup |5xi •••«9x„T(- -b (d-1) ^^‘^6-X))\ 
ee[o,i]" NeA’" r,t 


(3.3) 


where dy denotes the directional derivative in the direction Y, and for 6 £ [0,1]” and X £ 
we abbreviate 


e-x -.= 9iXi + --- + enXn. 














Proposition 3.1. Let H(t) be the constrained Dyson Brownian motion from Definition \2.S\ with 
initial condition (j2.4|) . Fix e > 0 and let r = r(e) he large enough depending on e. Then for any 
F £ C^{M) we have 


KF{H{t)) -KF{H{<d)) 


O 




max [ ||(9*F|Ls 



(3.4) 


In the applications in Section [U we will use functions F satisfying ^ for i ^ 4 

and a constant c > 0 that can be chosen arbitrarily small. Thus, for t ^ the 

right-hand side of (j3.4l) will be which is o(l) provided that c + e < 5. 

The starting point for the proof of Proposition 13.11 is the idea of |in( Lemma A.l], namely 
to estimate the left-hand side of (13.41) by estimating K{LF{H(t))). However, since the entries 
of H{t) are not independent, a different approach from [TO] is needed to control E(LF{F[{t))). 
We do this by approximating the constrained DBM by a Markovian jump process induced by 
switchings. This process is defined as follows. 

3.1. Switching dynamics. We introduce a Markovian jump process on simple regular graphs by 
defining its generator 


QfiA) ^ E (3-5) 

where we recall the definition of a switching from (12.5p and introduce the indicator function 

:= AjAmnil - Ai^){l - Ain){l - Aj^)il - . (3.6) 

The indicator function Ift^{A) ensures that the graph encoded by A contains the edges {i,j} and 
{m,n} but no other edges between the four vertices {i,j,m,n} (i.e. its restriction to {i,j,m,n} 
is 1-regular). 

Thus, the process generated by Q is a Markovian jump process whose jump times are the 
events of a Poisson clock with a constant rate; at each event of the clock, four vertices are 
selected uniformly at random, and a switching as in Figure 12.11 is performed on the graph if 
the four vertices are connected by exactly two edges. It is not hard to show that the uniform 
measure on d-regular graphs is invariant under this jump process. 

Proposition 3.2. The uniform measure on simple d-regular graphs is invariant under Q. This 
means that for any function f on the set of simple d-regular graphs we have K{Qf{A)) = 0. 

The proof of the proposition is given in Section [3.21 in a slightly more general context. The 
following proposition shows that the switching jump process generated by Q is well approximated 
by the constrained DBM generated by L. 

The generator L acts naturally on functions of H (denoted henceforth by an uppercase F), 
and the generator Q on functions of A (denoted henceforth by a lowercase /). It is therefore 
convenient to introduce, for any F £ C^(A4), the abbreviations 

H = Ha ■■= ^^(A-dee*), f{A) = fp{A) := F{Ha) ■ (3.7) 
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(3.8) 


Proposition 3.3. For any F G C^{A4) and using the notation (13.71) we have 


Qf{A) = LF{H) + R, 


where 

KR = ^a^^\d^F\\rfl . (3.9) 

Here E denotes expectation with respeet to the uniform measure on random d-regular graphs A. 

The proof of this proposition is also deferred to Section 13.21 below. Roughly, the idea of the 
proof is as follows. By Taylor expansion, we obtain 

Qf{A) » ^ E A,A„.„(-8.r'“/(A) + t(8’””)V(.4)) (3.10) 

with high probability. Now EAij = ^ if i / j. By expanding AijAmn = {jj + {Aj - w))(w + 
(Amn ~ ]v )), and keeping only the leading terms, we find that the right-hand side of (I3.10p 
becomes by LF{H). Here, for the second-order term on the right-hand side of (I3.1UI) . the 
leading term from AijAmn is for the first-order term on the right-hand side of (13.101) . the 
leading term from AijAmn is “ ;^) + ^{^mn ~ Further error terms result from the 

dependence of the entries of the adjacency matrix. 

Before giving the proofs of Propositions 13.2113.31 we deduce Proposition 13.11 from them. 

Proof of Proposition [37T1 By ()2.10p . it suffices to estimate E[LF{H{t))]. By explicit solu¬ 
tion of the constrained DBM, H{t), we find for any fixed t ^ 0 that 

H{t) = e"*/277(0) + (1-e"*)^/2lF (3.11) 


where IT is a copy of the constrained GOE independent of H{0). For the remainder of the 
proof, we identify the right-hand side with H{t), abbreviate H = H{0), and introduce the two 
functions 

Fw{H) = Fh{W) := F(e-*/^H +(1-e-*)^/2^) , 

where the choice of the argument determines the variables on which the generator L acts. We 
recall the generator L from (12.9p . 


L 


1 

16A^3 


E 'Mif 


^ \ ^ uklakl 

32iV2 ■ 

z,j,k,l 


From (9^ = (e * -|- (1 — e *))5^, e = dFw, and (1 — e = dFn, we then deduce 

that LF{e-^^^H + (1 - e-y/'^W) = LFw{H) + LFh{W). We therefore get 

E[LF{H{t))] = E[LFw{H)]+E[LFh{W)] = E[LFw{H)], 


where in the second step we used that the constrained GOE, IT, is invariant with respect to the 
generator L. 

Next, we define fw^) ■= Fw{H) where H = Ha is defined as (13.71) . By Proposition 13.21 
the random d-regular graph A is invariant with respect to the generator Q, and Proposition 13.31 
therefore yields 

E[LFw{H)] = E[Qfw{A)] + 0{D-^/'^N^+^)ui&^\yFw\\rfl = ui&^\yF\y ■ 
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Figure 3.1. Given four vertices i^j, k, I with two edges between them, there are two possible switchings. 
By equipping the edges with directions, one of these two switchings can be selected canonically. 


Thus, with (I2.10I1 . we have shown that 

^nF{Hm = 0(Zl-V2ivi+-)max||a*F||,,i, 
dt 

and the claim follows by integrating over t. □ 

3.2. Proofs of Propositions [3.2113.31 Propositions I3.2fl3.3l concern switchings of regular graphs. 
Switchings played an important role in the proof of the local semicircle law for random regular 
graphs [3]. Here we use simple switchings instead of the double switchings needed in [3]. 

Given two disjoint edges of a regular graph such that the graph has no other edges between 
the vertices incident to these two edges, there are two possible switchings; see Figure 13.11 To 
specify one of these two switchings, it is convenient to assign to each of the edges to be switched 
a direction; there is then a canonical choice between the two possible switchings. We write ij 
for the edge {i,j} directed from i to j. 

We consider sets S of two directed edges of the complete graph, which we write in the form 
S = {ij,kl}. We denote by [S'] = {i,j,k,l} the set of vertices incident to the edges of S. For 
two such sets S and S', we define the indicator functions 


I{S) = I{S;A) ■■= 1(|[S]| = 4 and £'|[ 5 ] is 1-regnlar), (3.12) 

J(S,S') = J{S,S';A) ■.= l([S]n[S'] = 0 ), (3.13) 

where E = E{A) := : Aij = 1} is the set of (undirected) edges of the graph encoded by 

A, and E\b ■= {e £ E ■. e C. B} is the restriction of the graph E to the subset of vertices B. The 
indicator functions are illustrated in Figure 13.21 Note that jmn ^ AijAmnI{{ij,mn}) (recall 

O). 

Remark 3.4. The definitions (j3.12p - (|3.13p are similar to those given in [3l Section 6], with the 
following differences. First, the cnrrent set S consists of two directed edges instead of the three 
undirected edges in [3]. Because of the directions contained in the current set S, it effectively 
incorporates the extra parameter s of [3l Section 6]. Second, the edges in S are edges of the 
complete graph, and we do not assume that they are contained in some regular graph A; we will 
ultimately define the switching associated with the set S to act trivially unless S is contained 
in the edges E of the given graph. 

For a set 5 = {ij, kl} of two directed edges, we define the switching 


Ts{A) 


'A-^fj AI{S) = l,Aij = l,Aki = l 
< A + C’:j AI{S) = l,Aik = l,Aji = l 
A otherwise, 


(3.14) 
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V_ J 


Figure 3.2. In the left diagram, I{S) = 0 since the restricted graph is not 1-regular. In the right 
diagram, J{S, S') = 0 since the two sets of vertices intersect. 


where we recall the definition of from (12. In words, Ts{A) switches the edges S if they are 
contained in A and are switchable in the sense that the switching results again in a d-regular 
graph. Moreover, for S', S' as above, we define 


TsA^) 


Ts^{Ts{A)) if J(S, S') = 1 
A otherwise. 


(3.15) 


In words, T 5 ^ 5 /(A) switches the edges in S and S' if they are contained in A and the two 
switchings do not interfere with each other. 

Lemma 3.5. For any fixed S, S' we have A = Ts{A) and A = Ts^s'i^)- 

Proof. It is easy to check that Ts{A) is a d-regular graph if and only if A is. Moreover, 
Ts{Ts{A)) = A, so Ts is a bijection on the set of d-regular graphs. Since the distribution of A is 
uniform, we obtain A = Ts{A). The second claim follows similarly from Ts^s'{Ts,S'{^)) = A. □ 

Now Proposition 13.21 follows easily. 

Proof of Proposition 13.21 For any /, we get 




nAmAjnI{{ij, mny,A + e^'')f{A + ^P")) 
Y, F.{Ai,AmnI{{i3, mn}; A)f{A - 4™")) 

Y. E{inA)f{A - £ 5 ”)). 


where the first and last steps follows from the definition of dij", the second step from Lemma [3. 5 1 
and the third step from the exchangability of i,j,m,n and using I{S]A) = I{S‘,Ts{A)). This 
concludes the proof. □ 


For the proof of Proposition 13.31 we shall need estimates on the moments of entries of the ad¬ 
jacency matrix, as well as estimates on such moments restricted to low-probability events where 
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the indicator functions (I3.12l) - (j3.13p are zero. These estimates are collected in the following 
sequence of lemmas. 

The following two lemmas show that moments of the entries of the adjacency matrix behave 
roughly like those of an Erdds-Renyi graph. 

Lemma 3.6. Let b N and Then for any p € [l,iV] and q € 

[l,iVl \{ii,ji,...,ib,jb}, we have 

''' ^ibjb^pg) ~ > (3.16) 

where we use the convention • • • Ai^j^) = 1 if b = 0. 

Proof. Set X := Ai^j^ ■ ■ ■ and I := {h,ji ,... ,ib,jb,p}- Then, since = d for any 

p, we find for any q ^ I that 

1E(X) = i + 

n n^I n£l 

n^I 

where in the third step we used that XApn ^ 0 and in the last step that the law of A is 
invariant under permutation of vertices. Using that |/| ^ N/2 by assumption on b, the claim 


now follows. □ 

As a consequence of Lemma 13.61 we obtain the following explicit bounds. 

Lemma 3.7. Suppose that |{i, j, m, n}| = 4 — a and \{i,j, k, I, m,n,p,q}\ = 8 — b. Then 

/ r] \ 2 -1 a/2j 

E(AjAmn) = 0[-) , (3.17) 

E{AijAmnAkiApg) = 0(^—j . (3.18) 

Proof. Since = 0 for all s, we can assume that i ^ j, m ^ n, k ^ I, and p ^ q, and thus 
a ^ 2 and 6^4. Then (I3.17l) - (j3.18p follow easily from Lemma [T6l □ 


In the next two lemmas, we estimate moments restricted to low-probability events where the 
indicator functions (j3.12p - (l3.13jl vanish, i.e. we estimate the contribution of graphs A that are 
not switchable. Throughout the rest of this section, for given indices z, j. A:,/, m, n,p, g we use 
the abbreviations 


h ■= I{{ij,mn};A), h := I{{kl,pq]\A), (3.19) 

Ji 2 ■= J{{ij,mn},{kl,pq}-,A), I 12 := hhJu, (3.20) 

with I and J defined in (j3.12p - ()3.13p . 

Lemma 3.8. Let |{z, j, m, n}| = 4 — a and |{i, j. A:,/, m, n,p, g}! = 8 — 6. Then 

^{{AijAmn + AimAjn){^ — II)) = ' ( 3 - 21 ) 

¥^{{AijAran AimAjn){AklApq -\-AkpAiq){\ — I 12 )) = ^ ' (3.22) 
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Proof. First, assume that i,j, k, I, m, n,p, q are all distinct, i.e. we consider the case a = b = 0. 
Then, since \{i,j,m,n}\ = 4 and /i = 0 implies that the graph A restricted to {i,j,m,n} is not 
1-regular, we hnd 


K{AijAmn{^ — Ii)) ^ K{AijAmn “1“ ^in AjYn “h , 

'^{AimAjn{l — II)) ^ 'K(AimAjn {Aij + Afyin T Ain T ; 

and Lemma [3.61 implies that the right-hand sides are bounded by 0{d/N)^. The proof of (I3.22p 
for 6 = 0 is analogous. We only consider the term AijAkiAmn-^pq] the others dealt with similarly. 
First, note that J 12 = 1 if |{i, j, k, I, m, n,p, (;}| =8. Since |{i, j, k, I, m, n,p,q}\ =8 and I 1 I 2 = 0 
imply that or has at least three edges, we hnd 


^(^{Aij AmnAf^lApq(^l Li/2>^12)) — 

^ K(^(AijAmnAklApq[Aim + Ain + + Ajn -|- Akp + A^q + Alp + Aiq)'j = > 

where the last step follows from Lemma 13.61 

Finally, if a > 0 we have Ii = 0, and if 6 > 0 we have /12 = 0. In these cases, we can directly 
apply (j3.17p and (j3.18p . respectively, and the claim follows since 2 — [a/2j ^3 — aifa>0 and 
4 - [b/2\ ^5-6if6>0. □ 

As a consequence of Lemma 13.81 we obtain the following averaged estimates. 

Lemma 3.9. If |{i, j}| = 2 — a and \{i,j, k, /}| = 4 — b, then 

-Jp, + AimAjn){'^ — h)) = J (3.23) 

m,n 

iiEE ^{{AijAmn + AimAjn){AkiApq AkpAiq){\. — / 12 )) — jy) ’ (3.24) 

m,n p,q 

Moreover, 

Jp ^{{AijAnin + AiniAjn){^ — II)) = —^ , (3.25) 

Jp ^{{AijAmn A AitnAjn){AklApq + AkpAiq)(l — I 12 )) = ' (3.26) 

i,j,m,n k,l,p,q 

Proof. To prove (13.231) . we split the summation over m,n by hxing \{i, j,m,n}\ = 4 — a — s 
where s £ [0,2]; there are 0{N‘^~^) terms corresponding to each s G [0,2]. By (I3.2ip . the 
left-hand side of (|3.23p is bounded by 


O 


d \3-a 


N 


+ 


Y,0iN-^)0[ 


S = 1 


d\3-a-s 


N 


= O 


d \3-“ 

N. 


The proofs of (I3.24ll - (j3.26jl are analogous. 


□ 
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Finally, as a consequence of Lemmas 13.6H3.9I and the Holder inequality, we obtain the follow¬ 
ing estimates incorporating an arbitrary function f{A). These and the remainder of the proof of 
Proposition 13.31 are simplest to state in terms of versions of the seminorms (I3.2p - (l3.3p for t = 0 
without rescaling by (d — Thus, instead of (j3.2p and (j3.3p . we use the seminorms 

ll/ll, := (E|/(/l)|’')'''’' 


and 


WdVWr 


sup sup \dxi ■■■dxr,f{ - + 0 - X)\ 

9&[0,1Y 


Lemma 3.10. Fix e > 0 and let r = r(e) be large enough depending on e. Let f G C^{M.) satisfy 
\\f\\r ^ 1- Then if |{i, j}| = 2 — a and |{i, j, k,l}\ =4 — 6, we have 


j^2 y^.^i^ij^mnfjA)) 

m,n 

m,n^p^q 

m^n 

wEE ^{{AijAyyin + AijYLAjn){A}^iApq A}^pAiq)Ii2f{A)^ 
m,n p,q 

^{{AijAmn + AijnAjn)Ilf{A)^ 
^{fAijAjyin T AiYnAjn){AkiApq Ai«pAiq)Ii2f {A)^ 

i,j,m,n k,l,p,q 


o 

o 

o 

o 

o 

o 


^ 2-la/2\-e 

n) 

d n4-L6/2J-£ 

n) 

d \ 3-a-e 

n) 

d \ 5-fe-e 

n) 

d \ 3-e: 

n) 

d \ 5-e 

n) 


(3.27) 

(3.28) 

(3.29) 

(3.30) 

(3.31) 

(3.32) 


where Ii := 1 — I\, I 12 := 1 — I 12 , and the indicator functions Ii and I 12 were defined in 

(l3T9D ^ (iM]l . 

Proof. We only prove (I3.3ip : the other estimates are proved similarly and we comment on the 
differences at the end of the proof. By Holder’s inequality, applied twice, first to E(-) and then 
to the sum over m, n, we obtain from (I3.25P that 

^{{AijAmnFAimAjn){i—Il)f{A)^ 

^ W [E((H,,H^„ + H,^H,„)(l-/i))]'-'/’'||/llr- 



where we chose r large enough that 3/r ^ e. 

To prove (13.321) . we use (13.261) instead of (I3.25p . and to prove p3.27l) - (|3.28p we apply (13.161) 
instead of (j3.2ip . To prove (I3.29j) - (I3.30I) . we use (j3.23p - (l3.24p . This concludes the proof. □ 
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The next lemma estimates the effect of replacing Aij by its mean d/N, or, equivalently, of 
conditioning on {Aij = 1}. Since the entries of A are not independent, we use switchings to 
analyse such a conditioning. 

Lemma 3.11. Fix e > 0 and let r = r(e) be large enough depending on e. For any f G C‘^{Ai) 
and any i,j, k, I with |{i, j}| = 2 — a and |{i, j, k,l}\ =4 — 6, we have 


E 



O 

O 


d \ i-£ 

n) 

_d\2-£ 

n) 


II5/II. + 0 

\\d^f\\r + 0 


d \2-a-e 

n) 

d \3-b-s 

N 


(3.33) 

(3.34) 


Proof. We begin with (I3.33p . Since A ^ A4 + dee*, we have Ylm n ^mn = Fid and Aim = 
Ajn = d for all i, and the left-hand side of ()3.33p is therefore equal to 


e(/(A)(a^-4)) = ^^e(/(/1)(A.,A 

\ \ / / lyyi 'n 


mn -^imAjn)'j 


(3.35) 


Using (I3.29p . using the notation from (I3.19D . we therefore find 

E(fiA)U,-^']'] = ±^E{fiA)iA,jAmn-A,mA,n)h)+0[^y~''~"\\f\\r. 

^ ^ ^ ^ m,n 


Because of the indicator function Ii, the first term on the right-hand side vanishes unless a = 0. 
Therefore we may assume that a = 0 when estimating it. By Lemma 13.51 and since Ii{A) = 
Ii{Ts{A)) with S = {ij,mn}, the first term on the right-hand side equals 

^ ^e((/(^) - fiA - ^^-))A,AmnIl) . (3.36) 

m^n 


The difference of the /’s is bounded in absolute value by supggjQ sup^^-g^ri^ \dxf{A+6X)\. Hence, 
(j3.27p implies that (|3.36p is bounded by 

o(v) 

This concludes the proof of (13.331) . 

The proof of (I3.34p is similar. As in (I3.35p . we write 




1 

{Ndf 


E 

m,n,p,q 


{Aij A. 


AimAjn){Al,,lA 


pq 


AkpAiq') . 


As above, we write 1 = I 12 + (1 — / 12 ) inside the expectation on the left-hand side of (I3.34p . 
The second term yields a contribution of order 0{-^)^~^~^\\f\\r, by (I3.30p . The first term is 
zero unless 6 = 0 because of the factor J 12 in I 12 . We may therefore assume that 6 = 0 for the 
estimate of the first term. Using Lemma 13.51 as in (I3.36p . we find that the first term is equal to 

^ ' m,n,p,q 

(3.37) 
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The difference of the four /’s is bounded in absolute value by 

sup sup \dx:,dx 2 f{A + 9iXi + 62 X 2)1 . 

By (j.S.28jl . we therefore find that (j3..S7p is bounded in absolute value by 



This concludes the proof. □ 


Finally, with the preparations provided by the previous lemmas, we now complete the proof 
of Proposition 13.31 


Proof of Proposition 13.31 First note that = AijAmnh- By Taylor expansion, and 
writing Ii = 1 + (/i — 1), we therefore have 

Qf{A] = ^ E + R,) , (3.38) 


where 


= - h) sup sup \dxf{A + 9X)\ 


Ri = O 




ee[o,i]Xex 


= 


R 2 = o 


dJ m 


X AijAmn sup sup |(9xi/(-4 + 6» • X) | . 

* ^ — itn t.r_n>Q 




6»e[o,i]3 xeT3 


By (I3.3ip and ()3.27p . respectively, the two error terms are estimated by 


Next, we estimate the main terms in (j3.38l) . which we write as 

1 


8 Nd 


A.,Au{-(f‘f{A) + l{d»ff(A))) 




(3.39) 


The idea is to write Aij = ^ + [Aij — and likewise for Aki- For the second-order term 
in (j3.39p . the term obtained by selecting both factors ^ yields the main contribution. More 
precisely, we write 


where 


^ a„a,M‘?j(a) = ^ J^(a«)y(yi) + jv^{7;3 + i;,) 


R 3 


i?4 


N 1 


i,j,k,l 




d \ d 


N 1 


i,j,k,l 


E(((Sgt/(.4)) Ay 


N J N 
N 


Aki 
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By (I3.33P and (|3.34l) . respectively, with / replaced by we obtain 

E{R, + R,) = 0[^y~"[\\d^f\\r + \\d^f\\r)+0[^y~"\\d^f\\r^ 

Next, we estimate the first-order term in (j3.39l) using a similar argument. Here the term 
obtained by selecting both factors ^ from Aij and A^i vanishes because ~ 

main contribution is given by the mixed term. More precisely, we write 


SNd 


Y, A,Akidf^f{A) 


" E + 7^ E ( 




where 


8 N^ 


i,j,k,l 


4iV2 




_ A 


d^Jf{A) + N^R5 


i,j,k,l 


E , 

i,j,k,l 


N 1 




= MjvjE((at'/('^))('^.i-v)(a'“-v 


N 


N 


By (j3.34p . with / replaced by d^jf, we obtain 


ER, = o[^y "\\d^f\\,+o[^y "wdfWr. 


We conclude that 


Qf{A) 


d 

16iV3 




d- 1 

16iV3 




Vd^ 

~32W 

Vd^ 

~32W 


2=1 

6 

ij,k,l 2=1 


where we defined 


Rq 


1 1 

imW 




Clearly, EHs = 0(^)||aVllr. 

Using the notations introduced in (|3.7p . we have y/d — 1 df{A) = dF{H). Hence we obtain 
(13.81) with R := N"^ Ri- The error term R is estimated, using the above estimates on ERi, 
as 


m = o(jv2+') (4)"(||8/||, + i|9VW + 4||aVllr + 4(i|ayi|, + ||a"/y 


= 0(E)-^/2jy1+£) 


||9H||,,o + D-^/^\\d‘^F\\r,o + ||53f||,,o + D-^/^\\d^F\\r,o 


as claimed. 


□ 
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4. Stability of eigenvectors and eigenvalnes 

In this section we derive basic stability properties for the eigenvalues and eigenvectors of the 
Dyson Brownian motion H{t). These allow us to deduce estimates on the eigenvalues and 
eigenvectors of H{t), assuming similar estimates have been proved for H{0). 

As discussed in Section EH we consider a general Dyson Brownian motion H{t) on an M- 
dimensional Hilbert space V, with normalization constant N as in (I2.1j) . For the usual DBM 
we have N = M, while for the constrained DBM we have M = N — 1; we always assume that 
N and M are comparable. We denote by Ai(t) ^ ^ ^M{t) the eigenvalues of H{t), and by 

vi(t),.. . S V the associated normalized eigenvectors of Moreover, we define the 

Stieltjes transform of the empirical spectral measure of ff(t) by s(t; z) := ^ x (t)-z • 

Throughout the rest of the paper, we use the following notion of high probability events and 
high probability bounds, introduced in HZ]. 

Definition 4.1. fij We say that an event H has high probability if for every C > 0 there 
exists an A^o(C) > 0 such that P(H'^) ^ N~‘^ for N ^ A^o(C)- 

(a) For nonnegative random variables A,B, we write A < B or A = 0^{B) if for any C > 0 
there exists an Nq{Q) such that P(A > ^ N~‘^ for N ^ No{(). 

If the event H from (i) and the random variables A and B from (ii) depend on some additional 
parameter u ^ U in some possibly N-dependent set U, we we say that (i) and (ii) hold uniformly 
in u if Nq{() does not depend on u. 

Throughout the following, the definitions (i) and (ii) will always be uniform in all parameters, 
such as z, any matrix indices, and deterministic vectors. Note that -< is compatible with the 
usual algebraic operations, so that for instance we have Hi ^ Si provided that Ai -< Bi 
for all i and the size of the index set for i is . 

4.1. Delocalization of eigenvectors. The following result shows that if all eigenvectors of H{0) 
are uniformly delocalized in some direction q G H, then with high probability they remain 
delocalized in this direction under the DBM on V, for any time t > 0. 

Lemma 4.2. Suppose that H{t) is the DBM on an M-dimensional space V. Let q G H and 
suppose that maxj |q • Vj(0)| ^ B. Then, for any t > 0, any i G |1, M], and ^ 1, 

P(|q-Vi(f)| > ^ (4.1) 

In particular, 

|q-Vj(t)| -< B. (4.2) 

Lemma [4. 2 1 is a simple consequence of the eigenvector moment flow (EMF) introduced in [lOj . 
Suppose for simplicity that the eigenvalues of F1(0) are distinct. Then the eigenvalue process 
(Aj(t)) is almost surely continuous and simple for all t > 0; see m for more details. We study 
the dynamics of the eigenvectors Vj(t) by conditioning on the eigenvalue process; see again [TO] 
for a precise construction. Hence, for the following argument, we condition on (Aj(t)) and regard 
the eigenvalue process as deterministic. 

We give the definition of the EMF restricted to moments of a fixed order p G N. The 
configuration space is Clp := {p = {r]i))fi G : H^i =p}- The conhgurations rj £ Lip 
are interpreted as configurations of p particles on the lattice |1,M], whereby a single site of 
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[1, M] may be occupied by multiple particles. We denote by 77 *’-^ := rj + lirji > 0)(ej — ej) the 
configuration obtained from rj by moving one particle from i to j. The time-dependent generator 
R{t) of the EMF is defined by 

{R{t)f){r]) := - f{r ])), 

where ^ 

iV(A,(t)-A,(t))2' 

For our purposes, the precise form of the coefficients Wij{t) is not important; we only use 
that they are nonnegative and continuous in t. The p-particle FMF is given by the equation 


dtftiv) = {R{t)ft){r ]), fo:Qp^R given. (4.3) 


This is a linear (time-dependent) ODF on a finite dimensional vector space, and thus well-posed. 
It is also easy to see that it is contractive on in the sense that \\ft\\L°°(np) ^ ll/o||L°°(r 2 p)- 

Next, for deterministic rj € and q G V, we define 




:= E 


■ M 

n 

.i=l 


1 


(2r?i - 1)!! 


(q-v*(t))^^* 


(Ai(t) : i G |l,M],t ^ 0 ) , 


(4.4) 


where n!! := n • (n — 2) • • • 3 • 1 for odd n, and by convention (—1)!! = 1. In [lOl Theorem 3.1] it 
is shown that ft solves (|4.3I1 . 


Remark 4.3. In m, Dyson Brownian motion is defined without the Ornstein-Uhlenbeck drift 
term in the SDF (|2.2p . and the SDFs for the eigenvalues and eigenvectors are stated in m 
Definition 2.2]. In the present case, with drift term, the SDFs for eigenvalue and eigenvector 
flows are given by 


dAi 

dvj 




dBji 

y/N 

1 


y /. A^ Xj 




dR 




1 




dt 


2N (Ai - XjY 


for i = 1,2,..., M, and with B{t) a Brownian motion on the space of M x M real symmetric 
matrices with quadratic covariation {Bij, Bki){t) = {6ikSji + Sii6jk)t. Thus, the SDFs for the 
eigenvectors are the same with or without the drift term. Therefore the arguments of m 
Section 3] apply verbatim in our setting as well. 

Proof of Lemma 14.21 Suppose first that H (0) has simple spectrum. Let ft be the given by 
(I13D, which solves (|4.3p as remarked above. Then, since the EMF (14.3p is a contraction on 
L°°(Dp), we obtain from the assnmption of Lemma 14.21 that 

max 1 / 4 ( 77 )] ^ max 1 / 0 ( 77)1 ^ B'^p . 


Therefore, choosing rj = pet, we get 

E[(q.V4(t))2P] = (2p-l)!!E[/4(77)] ^ {2p - 1)\\B^P , 

from which the claim follows. Finally, if H{0) does not have simple spectrum, the same estimate 
holds by a simple approximation argument using the continuity of the eigenvectors as functions 
of the matrix. □ 
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4.2. Stability of eigenvalues. The following result shows that if the empirical spectral measure 
at t = 0 is close to the semicircle law, this remains true for t > 0. For its statement, recall that 
s{t, z) denotes the Stieltjes transform of the empirical spectral measure of H{t). We denote the 
Stieltjes transform of the semicircle law by m. It can be characterized as the unique holomorphic 
function m : C+ —)• C+ such that m? + mz + 1 = 0 and m{z) ^ 1/z as \z\ ^ oo; see e.g. [2]. 

Lemma 4.4. Suppose that C~^M ^ N ^ CM. Fix e > 0. If for some B ^ N ~^ we have 


s(0; z) — m{z)\ + i? + 


1 

(A^r/)V4 


(4.5) 


uniformly for z = E + ir] with p ^ N , then for any t ^ B we have 


s{t;z) 


m{z)\ + i? + 


1 

(lVr?)V4 • 


(4.6) 


uniformly for z = E + ir] with rj £ [N , 1]. 


Proof. Define Sfc^tiz) as the unique solution C+ —)• C+ of the self-consistent equation 


SfcA^) 


^ e */2Ai(0) - z - (1 - e ^)sic^t{z) 


(4.7) 


Thus, Sfc,t(z) is the Stieltjes transform of the free convolution of the empirical eigenvalue dis¬ 
tribution of e“*/^L7(0) and the semicircle law rescaled by (1 — e“*)^/^. We refer to [5] for the 
existence and uniqueness of Sfc,t( 2 :) and relative properties on the free convolution with semicircle 
law. 

As in (jd.lljl . we find that H{t) = e“*/^Lf(0) + (1 — e“*)^/^lT, where W is the standard 
Gaussian measure on T-L{V) with inner product (12.ip . Under the assumptions of the lemma, |33l 
Corollary 7.11] implies that for t ^ N~^ we have 

|s(f;z)-sfc,t(2)| + (4.8) 

uniformly for z = E + 17] with p ^ (Note that in [33], the Stieltjes transform is denoted 

by mv instead of s, and that is denoted m^ct- Moreover, m Corollary 7.11] is stated for 
a diagonal matrix H{0); however, since W is invariant under orthogonal transformations which 
diagonalize 77(0), the results of [33| trivially apply to any symmetric matrix 77(0).) 

Set Dt •= 1—^ t. Note that the Stieltjes transform of the empirical eigenvalue distribution 
of e“*/^77(0) is given by e*/^s(0, e*/^z), and that (14.7p can be rephrased as 

Sfc,t(z) = e*/^s(0,e‘/^(2+ i?tSfc,t(2;))) . 

For any z = E + ir] such that t] ^ we have Ime*/^(z + 'dtSfc,t{z)) ^ Ime^/^z ^ , 

where we used that Imsfc^t(z) > 0. From the assumption ()4.5p we therefore get 

Sic,t{z) = e^/'^m{e^/'^{z + ^tSic,t{z))) + • ( 4 - 9 ) 

Next, note that 

m{z) = e*/^m(e*/^( 2 : + dtm( 2 ;))). (4-10) 
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(This may be interpreted as the fact that the semicircle law with variance t is a semigroup with 
respect to free convolution.) Moreover, from the dehnition of m{z) it is easy to deduce the 
continuity estimate 


m{z) 


Tn{w)\ ^ 2 \z — , 


(4.11) 


(A^r/)V4 


for any z,w £ C+. 

By (14.lip , and using that t = 0(1), the difference between ()4.9p and (|4.10p is 
- m{z)\ = e*/^ |m(e*/^(z +-dfSfc,t( 2 :))) - m(e*/^(z + ??tm(z))) +0^(^B + 

^ - m(z)|^/" + 0^(^B + 

^ max|o(t^/2^|sfc,t(z)+ . 

Therefore either |sfc,t( 2 ) — m(z)\ = 0{t) or |sfc,t(z) — m(z)\ -< B + and we get 


\sfc,ti^) - m{z)\ ~< B+ 


+ t ^ B + 


(4.12) 


(A^r/)V4 (Afr/)V4 ’ 

where we used t ^ B. Combining (14.8p and (|4.12p and using rj ^ 1, the claim (j4.6p follows. □ 


5. Proof of Propositions 12.61 - 12.71 

With the preparations provided by Sections [SHU and using results of [311271133] , we now complete 
the proofs of Propositions 12.6112.71 First, recall that a > 0 is fixed, and that we always assume 
D ^ N°‘. We also use the notation z = E + \r] for the real and imaginary parts of the spectral 
parameter 2 : S C+. 

Throughout this section, H{t) denotes the constrained DBM from Definition 12.21 with 77(0) 
given by (|2.4p . We use the notations of Section [4] applied to the constrained DBM. In particular, 

M := N-1 


is the dimension of the space V := e-*-. 


5.1. A priori estimates on eigenvalues and eigenvectors. We begin by collecting some results 
on the eigenvalues and eigenvectors of H{t)\^± required for the proofs of Propositions I2.6H2.71 
For any H G A4, we denote the eigenvalues of i7|e-L by Xi{H) ^ ^ Xm{H), and the corre¬ 

sponding normalized eigenvectors by vi(i7),..., vm{H). The components of the eigenvectors in 
the standard basis on are denoted Vk{H]i) := e* • Vfc(77), i € |1, A^J, k G [1,M]. Moreover, 
for H G A4, we denote by Gij{H\z) the entries of the Green’s function of H restricted to e-*- 
in the standard basis of and by s{H;z) the Stieltjes transform of the empirical spectral 
measure. Explicitly, 


Gi,{H-z) 


^ Vk{H]i)vk{H;j) 


s{H]z) 


M 


TtG{H;z) 


1 

M 


E 


1 

XkiH) - z ■ 


(5.1) 

(5.2) 
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Finally, we set 


(5.3) 


r{H) = r{H-,z) := max\Gii{H;z)\\/l. 

We also recall the definition of the typical location 7 ./ of the i-th eigenvalue from (II.ip . 

The following proposition summarizes the input we need from the local semicircle law of [3]. 
The local semicircle law, as proved in [3], only applies for t = 0, and the extension to t > 0 is 
provided by the results of Section 01 

Proposition 5.1. For any z G C+, i G |1, k G and 0 ^ ^ we have 

\vk{H{ty,i)\ X nH{ty,z) ^1 + ^- (5.4) 

Moreover, for any fixed k > 0 and any i G (1 — we also have 

(5.5) 

Proof. First, as special cases of [H Theorem 1.1 and Corollary 1.2], for any z = E + \r] with 
FI G M and rj for arbitrary e > 0, we have 

\s{H{0yz)-m{z)\ < ^ + \vk{H{Qyi)\ < N-^/\ (5.6) 

(Note that the local semicircle law from [3] also includes the trivial eigenvalue at 0; it is easy to 
see that its contribution to s is negligible compared to the error bounds in (I5.6l) .l 

Next, we extend these bounds from t = 0 to t > 0. For i G |1, A^J define e* = e* —(ej-e)e G e-*-. 
Since Vk{H{t)-,i) = e* • Vfc(Ff(t)), from (|5.6I) and Lemma 14.21 applied to the constrained DBM 
with q = Oj, we find \vk{H (t)] i)\ -< for any t > 0. Similarly, for t ^ the extension 

of the bound on the Stieltjes transform follows immediately from Lemma 14.41 with B = 
Summarizing, for any rj ^ and 0 ^ ^ we have 

\s{H{tyz) -m{z)\ < \vk{H{tyi)\ < . (5.7) 

This proves the first estimate of (15.41) . 

In order to prove the second estimate of (|5.4p . we use a dyadic decomposition (see e.g. [25l 
(8.2)]) to obtain, for any matrix H G Ad, 

M . / riog2»7"P 

\Gij{z)\ ^ \Xk-^E + irj\ ^ 1 + ^ Ims(.F + i2”7/) 

We apply this estimate to the matrix H{t). By (15.71) . we have max^^/ \vk{l)\'^ -< l/N. Moreover, 
since rjlms{E + iij) is increasing in rj (as may be easily seen from the right-hand side of (15.21) 1. 
and since jmj ^ 1, the first bound in (|5.7p implies Ims( 2 ;) -< 1 -|- l/{Nrj) for any rj > 0, and thus 
Ims(Fl -|- i2"r/) -< 1 -|- 2~'^ /Nrj. For rj we then have log 77 “^ ^ 1 and obtain F(z) ^ 1 

as desired. For arbitrary rj > 0 the claim then follows by [3l Lemma 2.1]. (In fact, we shall only 
need (|5.4I) with rj 

Finally, we deduce (15.5|) from the bound on the Stieltjes transform in (15.71) . We abbreviate 
Afc = Xk{H{t)), and denote by 

Qik) ■= / g{x)dx, v{I) := ^ I) 

^ U — -1 
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the semicircle and empirical spectral measures, respectively, applied to an interval I. Then, 
following a standard application of the Helffer-Sjostrand functional calculus along the lines 
of |19l Section 8.1], we find from ()5.7p and D ^ N that for any interval I C [—3,3] we have 

W{I)-g{I)\ -< + -< (5.8) 

(We note that previously (|5.8p for t = 0 was given in [3], Corollary 1.3].) Using (15.81) . we may 
estimate A* — 7* as follows. By (15.81) applied to / = [—3,3], we find that there are at most 
eigenvalues outside [—3,3]. Defining f{E) := ^([£',00)), we therefore hnd from 
(HH) and (15. 8p that 

fin) = ^ = ^([^ooo))+ o(^) = z/([Ai,3))+ 

= e([Ai,3))= /(A*)+ 0^(-^) . 

Since i G |kA^, (1 — ^iV)], we have [/'[ ^ c > 0 in a neighbourhood of 7^, and we therefore get 
(|5.5I) . This concludes the proof. □ 

The next result shows that the suprema in (|3.3I) do not essentially change the size of T. 

Corollary 5.2. Fix n G N. For any z G C+ and 0 ^ t ^ , we have 

sup sup T(H{t) + {d—l)~^^‘^9-X-,z) -< 1 + -^. (5.9) 

0e[o,i]"XeY" Nr] 

Moreover, for any i G |1, W] and k G |1,M], we have 

sup sup \vkiH{t) F {d — 1)~^N ff . x-,i)\ f _ (5.10) 

6»e[o,i]'* xeY" 

Proof. We abbreviate H = H{t). Without loss of generality, by an argument analogous 
to [U Lemma 2.1], we may assume that rj ^ 1/N. Hence, by (15.4p . we have T{F[]z) -< 1. It 
therefore suffices to show that if T{F[] z) ^ {d— l)^/^/(16n) then for any 6 G [0, !]"■ and X G T” 
we have 

T{H+{d-l)-^N0.X;z) F 2T{H-z). (5.11) 

To show (15.lip , we use the resolvent identity to obtain (omitting the argument z for brevity) 

\Gi,{H + {d-l)-^/^e-X)\ 

= -id- l)-^N(^GiH)i6 . X)G{H + (d - 1 )-^N 6 ■ X)^ ^ 

^ T{H) + 8n(d - l)-^NriH)r{H + (d - l)-^/^ q . x) 

F ViH) + r(£ + {d- l)"^/2 e-X)l2 . 

Taking the maximum over i and 7 yields (15.111) . Finally, (15.101) follows from (15.91) . as in the 
proof of la Corollary 1.2]. □ 

Note that since GijiH;z) = GijiH;z), the estimates (15. 4p and (15.91) for T also hold with 
77 < 0 if 77 is replaced by I77I on the right-hand sides. We shall use this tacitly in the following. 
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5.2. Proof of Proposition 12.6t eigenvalue correlation functions. We now prove that the locally 
averaged local correlation functions of the matrix ff(0)|g± converge to those of for times 

t ^ . The main ingredient of the proof is the following lemma comparing functions of 

Green’s functions with spectral parameter rj slightly smaller than 1/N. Its proof follows easily 
from Proposition 13.II and Lemma [5.II For random matrices with independent entries, analogous 
results were previously proved by the Green’s function comparison theorem [25], and by direct 
analysis of the evolution of the matrix entries under Dyson Brownian motion m- We also 
remark that, in [44|, eigenvalues are compared directly without involving the Green’s function. 

Lemma 5.3. Fix n G N, and let K,'y,6 > 0 be suffieiently small. Then the following holds 
for any r] G any sequence of positive integers ki,k 2 , ■ ■ ■ ,kn, any set of eomplex 

parameters zT = ± irj, where j G ^ ^ 2 — k, and the ± signs are 

arbitrary. Let (j) G he a test funetion such that for any multi-index m = (mi, • • • ,mn) 

with 1 ^ \m\ ^ 4 and for any a; > 0 sufficiently small, 

max{|9"^<^(x)| : \x\ ^ N^} ^ , (5.12) 

max||(9™’(/)(x)| : |x| ^ ^ . (5.13) 

Then, with the notations Gi{z) '.= z) andG 2 {z) ■= G{H{t); z), for any t ^ , 

we have 

Ef Tt Gi(z]) j ,..., AT-^" Tr Gi{z^)^ j - EfiGi ^ G 2 ) 

= 0(iV-‘5/2+®(^)). (5.14) 

Here, 4>{Gi —)■ G 2 ) is the expression obtained from the one to its left by replaeing Gi with G 2 . 
The implicit constants depend on n, ki,... ,kn, mi,...,m^, and the constants in (I5.12p - (I5.13I) . 

Proof. For simplicity of notation, we show ()5.14p only for n = 1 and ki = 1; the general case 
is analogous. We then write z instead of z\. To show the claim, it then suffices to show that 

lEf {N-^ Tv G{H(t); z)) -Ec/){N-^Tt G{H{Of z))\ = 0(tD-^/2jYi+5/2^o(7)) _ ( 515 ) 

Set F{H) := (l){N~^ Tv G{H-, z)). We claim that if r and n are fixed (arbitrarily, independently 
of N), and if t ^ for any sufficiently large N (depending on r,n,5), we have 

sup \\d'^F\\r,s ^ N^/^+oG) . (5.16) 

Given (15.161) . Proposition 13.11 with e = 6/4 yields (15.151) . 

Thus, it only remains to show (I5.16|) . Recall that the derivative of the Green’s function in 
the direction of a matrix X ^ X \s given by dxG = —GXG (using that elements in X act on 
e-*-). Therefore, by the Leibniz rule, for any A"i,..., X^ G X and any FI G A4, we have 

dx^---dx^G = i-ir GX,(i)G---GX,(„)G, 

<y&S„ 

where Sn is the set of permutations of n elements, and we omit the dependence on H on both 
sides in our notation. Since (with respect to the standard basis of M'^) each X ^ X has at most 
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8 nonvanishing entries, and since these are in {±1}, by definition of F it follows that 


N 


\N-^Tvdx,---dx^G\ ^ lV-iJ]n!max (GX,(i)G• • • 

CFGOti 


2=1 


i-pn+l 


< 8”n!r 


From this and the chain rule, we obtain that there exist constants Cn such that 

|5xi---9x„</>(iV-'TrG)| ^ C^F^- max 


(5.17) 


By Corollary I5.2l and since \r]\ ^ N ^ , we have supggjg supj 5 ('g;^;>n T{H{s) + {d—l) ^Gq-X) -< 

N"', for any 0 ^ s ^ t. For n ^ 4, by assumption (I5.12p and (15.171) therefore 

sup sup \dxx--dxr.(p{^~^^^G{H{s) + {d-l)-^/^9-X))\ ^ , ( 5 , 18 ) 

6»e[o,i]'" xeA’" 

On the complement of the high-probability event of ^ in ()5.18p . we use the trivial bound 
F ^ r]~^ ^ and (|5.13p . We obtain 

sup sup \dxx--dx„(l>{N-^TiG{H{s) + {d-l)-^G0.x))\ ^ ^ , 

6»e[o,i]'" xeA’" 

(5.19) 

for any 0 ^ s ^ t. By combining the estimates (I5.18h - (I5.19I) . for any constant r = 0(1), we 
have 

^ N^/i+OG) + ^ jV^/4+o(7)^ (5^20) 

where C is as in Definition 10 and chosen sufficiently large, depending on r. This concludes the 
proof. □ 


The following lemma is essentially |25p Theorem 6.4]. It transforms the statement about the 
Green’s function of Lemma 15.31 to a statement about the local correlation functions. 

Lemma 5.4. Consider two random matrix ensembles Hi and H 2 with Green’s functions Gi{z) 
and G 2 {z). Suppose that, for all cf and parameters as in the statement of Lemma\5f^ the estimate 
(I5.14P holds. Then the local bulk eigenvalue correlation functions of Hi and H 2 coincide. 

Proof of Proposition 12.6t correlation functions. The proof follows directly by com¬ 
bining Lemmas 15.3115.41 with 6 given as in the assumption of Proposition 12.61 □ 


5.3. Proof of Proposition [2T6l eigenvalue gap statistics. To prove that the eigenvalue gap statis¬ 
tics are stable for short times, we require a weak level repulsion estimate. Such an estimate was 
derived in 1271 Theorem 4.1] for sparse matrices with independent entries, using a level repulsion 
estimate for t G established in [33]. Here we adapt the proof of m Theorem 4.1] to 

random regular graphs. The nontrivial dependence is dealt with by Proposition 13.11 
If \i{H) is a simple eigenvalue of H\^±, we define 


Q^{H) 


1 

N2 


E 


1 

(A,(77)-A,(H))2’ 


(5.21) 


and extend this definition by Qi{H) := 00 if Xi{H) is not a simple eigenvalue. This quantity 
plays an important role in [33], where it is observed that it captures the singularities of the 
derivatives of A* (77). In [27], it is found that Qi is stable under DBM and can thus be used to 
show weak level repulsion from such an estimate for larger times (when a Gaussian component 
is present). 
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Proposition 5.5 (Level repulsion). Fix k > 0. Then for any sufficiently small r > 0, any 
i G {kN, (1 — k)A^|, and any s ^ 0, we have 

F{Q^{H{s)) ^ = 0{N-^/^). (5.22) 

In particular, 

F{\{H{s)) - \+i{H{s)) < iV-i-") = 0{N--/^). (5.23) 

Proof. The proof is analogous to that of 1271 Theorem 4.1], with H\^± instead of H. We here 
focus on the differences. These result from the replacement of |27l Lemma 4.3] by Proposition l3.ll 
which takes into account the nontrivial correlation structure of the random regular graph. As 
in m, if \i{H) is a simple eigenvalue of H\^±, we define the matrix 


Ri{H) 


E 


1 


Xi{H) - X,{H) 




^ r G{H-z) 

2m ^i{H) - z 


where oj is chosen such that the contour |z — Xi{H)\ = uj encloses only Xi{H). Then we have 

Q,{H) = 

Given r > 0, define a cutoff function x satisfying the following two properties: (1) y is 
smooth, and the first four derivatives are bounded, i.e. |x^^^(^)l = 0(1), for k = 1,2,3,4; (2) 
On the interval [0, |x(®) — 2;| ^ 1, and for x ^ x(^) = Then x° Qi extends to 

a smooth function on the space of symmetric matrices. 

The proof of (j5.22|l consists of three steps. The first step is the estimate 

E[x(g*(77(s)))] = 0(iV302) ^ (5.24) 

for s ^ t ■= This estimate follows from [331 Theorem 3.6], whose assumptions are 

satisfied with high probability for the random d-regular graph by Proposition 15. II In particular, 
independence of the entries of H is not used. 

In the second step, we derive the comparison estimate 


|E[x(g.(/7(t)))]-E[x(Q.(77(s)))]| ^ 1, (5.25) 

for s G [0,t]. Instead of using [271 Lemma 4.3], which requires that the entries of the random 
matrix H{s) are independent, we use Proposition 13.11 which takes into account the nontrivial 
correlation structure of the random regular graph. By Proposition 13.II with F{H) := x(Qi(77)), 
it suffices to bound 


Wd'^F 


r,s 


= E 


sup sup 
ee[o,i]" 


dx,dx, • • • dx,,F{His) + (d - 1)-^^ 9 • X) 



(5.26) 


for any (large) fixed integer r and n = 1, 2,3,4. To this end, the computation of the proof of [271 
Proposition 4.6] applies, by simply replacing the derivatives by dx^ ■ ■ ■ with Xi G A. 
Here the formulas [271 (4.16)-(4.18)] remain valid after replacing V by the Xi appropriately, 
and similarly the formula below [271 (4.18)] remains valid after replacing Vij by v*{H)XiVj{H). 
Moreover, an analogous formula holds for n = 4; see e.g. [391 p.8]. The same formulas are valid 
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with H replaced hy H + {d—1) 0 ■ X. Since the Xi have only 8 nonvanishing entries (in the 

standard basis on M'^), and these are equal to ±1, Corollary 15.21 then implies 

sup sup IV* {H{s) + {d- 9-X)Xi {H{s) + {d - 1)"^/^ 0 • X) | X iV"^ 

6»e[o,i]" 

for any s G [0,t]. As in the proof of \27\ Proposition 4.6], we therefore get 

sup sup \dx,dx^ ■ ■ ■ dx„F{H{s) + {d- 9 ■ X)\ ^ . 

6»e[o,i]'* x&x-^ 

From this, bounding (I5.26P as in (I5.20p . we obtain 

^ ^c+(n+2)r_ (5^27) 

for arbitrarily small c > 0 and N large enough. Then (|5.25l) follows from Proposition 13.41 since 
^ ^ 1 for t ^ and D ^ by chooosing c and r 

sufficiently small. 

In the last step, we combine (I5.24p and (I5.25p . and thus obtain 

E[x(Q*(77(s)))] = 0(X3-/2), 

for any s ^ 0. Then ()5.22p follows easily by Markov’s inequality and the definition of x- ^ 

Proof of Proposition [T6l gap statistics. Throughout the proof, we use the abbreviation 
\i{t) = Xi{H{t)). Fix K > 0, (5 > 0, and t ^ N~^~^D^/'^. Since ^( 7 i) is bounded above and 
below for i G |kX, (1 — «:)X], it suffices to prove (j2.14p with Q{yi) replaced by 1. Moreover, for 
any n G N and (j) ^ C°°(K"') with bounded first four derivatives, it suffices to show the stronger 
claim 

E(^(XAi(0),... ,XAi+„(0)) = E(j){NXi{t),... ,NXi+nit)) + o{l) (5.28) 

as X —>• oo, uniformly in i G {nN, (1 — k)X]. For simplicity of notation, we only prove (15.281) 
for n = 1; the general case is analogous and we comment on the differences at the end of the 
proof. Thus, for any i G [kX, (1 — k)NJ and (p G C'°°(M) with bounded first four derivatives, we 
show 

E[<))(XA(0))]-E[<))(XAi(t))] = o(l). (5.29) 

Given a small constant r > 0, we choose a cutoff function p such that p{x) = 1 for x ^ X^’’ 
and p{x) = 0 for x ^ 2X^’'. Using (|5.22p . we can first remove a bad event on which Qi is large: 

|E[<).(XA(O))]-E[0(XAi(t))]| 

^ \E[(P{NX,{0))p{Q,{H{0)))]-ncP{NXi{t))p{QiiHm]\ 

+ Halloo(nft(^(0)) ^ X2-) +P(Q,(X(t)) ^ X2-)) 

^ |E[().(XA(0))p(Q,(X(0)))] -E[<)>(XAi(t))p(Q,(X(t)))]| . 

To estimate the right-hand side, we apply Proposition 13.11 with F{H) := (j)(NXi{F[))p{Qi{F[)). 
By an argument analogous to that used to obtain (15.271) . for any r and n = 1, 2, 3,4, we find the 
bound 

^ iV^+o(r) (5_3Q) 
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for arbitrarily small c > 0 (and N sufficiently large). More precisely, by the product rule, the 
derivatives act either on (p{NXi) or poQ^. In the bound of any of these derivatives, by definition 
of p, we can assume that Qi ^ 2N‘^'^. Then the derivatives of p o Qi are bounded exactly as 
in the proof of Proposition 15.51 For the derivatives of by the chain rule and since (p is 

smooth, it suffices to bound the derivatives of the eigenvalues A*. This is again done similarly to 
the bounds on the derivatives of Qt. Indeed, the derivatives of the eigenvalues can be expressed 
in terms of the eigenvalues and eigenvectors as done in [271 (4.16)-(4.18)] (and with [391 p.8] for 
n = 4). The latter expressions are bounded using the delocalization of eigenvectors (|5.4|) . and 
using that 


E 


^ NQy\H{s)), Y. 


|Ai(s) - Aj(s)| 




\Xi{s) - Xj{s)\’^ 


^ N^QYiHis)), 


as in [271 (4.11)-(4.12)]. 

As a consequence of Proposition 13.11 and (I5.30p . with t ^ we finally obtain 


\E[p{NXi{0))piQiiH{0)))]-E[<P{NXi{t))p{Q,{Hit)))]\ = 0{N' 


c+0{t)—5 


), 


and (I5.29|) then follows by taking c and r small enough that c + 0(t) < S. 

In the general case of a test function (j){NXi ,..., NXi+n), we use the product of cutoff func¬ 
tions (p o Qi) • • • {p o Qi+n) instead of p o Qj, and proceed otherwise analogously. □ 

5.4. Proof of Proposition 12.71 

Proof of Propositions ETl Given the estimates (|5.4p - (|5.5p . the same argument as in [271 
Section 3] applies. □ 


Acknowledgements 

AK was partly supported by Swiss National Science Foundation grant 144662. HTY was partly 
supported by NSF grant DMS-1307444 and a Simons Investigator fellowship. We thank Ben 
Landon for helpful discussions, Peter Sarnak for informing us about |28ll36j . and Michael Aizen- 
man for bringing the works |3nil31j to our attention. 


References 

[1] O. Ajanki, L. Erdos, and T. Kruger. Local eigenvalue statistics for random matrices with 
general short range correlations, 2016. 

[2] G.W. Anderson, A. Guionnet, and O. Zeitouni. An introduction to random matrices^ vol¬ 
ume 118 of Cambridge Studies in Advanced Mathematics. Gambridge University Press, 
Cambridge, 2010. 

[3] R. Bauerschmidt, A. Knowles, and H.-T. Yau. Local semicircle law for random regular 
graphs, 2015. Preprint, arXiv: 1503.08702. 

[4] F. Bekerman, A. Figalli, and A. Guionnet. Transport maps for /3-matrix models and uni¬ 
versality. Comm. Math. Phys., 338(2):589-619, 2015. 


29 











[5] P. Biane. On the free convolution with a semi-circular distribution. Indiana Univ. Math. 
J., 46(3):705-718, 1997. 

[6] P. Bleher and A. Its. Semiclassical asymptotics of orthogonal polynomials, Riemann-Hilbert 
problem, and universality in the matrix model. Ann. of Math. (2), 150(l):185-266, 1999. 

[7] O. Bohigas, M.J. Giannoni, and C. Schmit. Characterization of chaotic quantum spectra 
and universality of level fluctuation laws. Phys. Rev. Lett., 52:1-4, Jan 1984. 

[8] P. Bourgade, L. Erdos, and H.-T. Yau. Universality of general /3-ensembles. Duke Math. 
J., 163(6):1127-1190, 2014. 

[9] P. Bourgade, L. Erdos, H.-T. Yau, and J. Yin. Eixed energy universality for generalized 
wigner matrices. Commun. Pure Appl. Math., 2015. 

[10] P. Bourgade and H.-T. Yau. The eigenvector moment flow and local quantum unique 
ergodicity. Comm. Math. Phys., pages 1-48, 2016. 

[11] Z. Che. Universality of random matrices with correlated entries, 2016. 

[12] N.A. Cook. On the singularity of adjacency matrices for random regular digraphs, 2014. 
Preprint, arXiv:1411.0243. 

[13] C. Cooper, M. Dyer, and C. Greenhill. Sampling regular graphs and a peer-to-peer network. 
Combin. Probab. Comput., 16(4):557-593, 2007. 

[14] P. Deift, T. Kriecherbauer, K.T.-R. McLaughlin, S. Venakides, and X. Zhou. Uniform 
asymptotics for polynomials orthogonal with respect to varying exponential weights and 
applications to universality questions in random matrix theory. Comm. Pure Appl. Math., 
52(11);1335-1425, 1999. 

[15] I. Dumitriu and S. Pal. Sparse regular random graphs: spectral density and eigenvectors. 
Ann. Probab., 40(5):2197-2235, 2012. 

[16] F.J. Dyson. A Brownian-motion model for the eigenvalues of a random matrix. J. Mathe¬ 
matical Phys., 3:1191-1198, 1962. 

[17] L. Erdos, A. Knowles, and H.-T. Yau. Averaging fluctuations in resolvents of random band 
matrices. Ann. Henri Poincare, 14(8): 1837-1926, 2013. 

[18] L. Erdos, A. Knowles, H.-T. Yau, and J. Yin. Spectral statistics of Erdos-R&yi Graphs 
H: Eigenvalue spacing and the extreme eigenvalues. Comm. Math. Phys., 314(3):587-640, 
2012. 

[19] L. Erdos, A. Knowles, H.-T. Yau, and J. Yin. Spectral statistics of Erdos-R&yi graphs I: 
Local semicircle law. Ann. Probab., 41(3B):2279-2375, 2013. 

[20] L. Erdos, S. Peche, J.A. Ramirez, B. Schlein, and H.-T. Yau. Bulk universality for Wigner 
matrices. Comm. Pure Appl. Math., 63(7):895-925, 2010. 

[21] L. Erdos, J.A. Ramirez, B. Schlein, and H.-T. Yau. Universality of sine-kernel for Wigner 
matrices with a small Gaussian perturbation. Electron. J. Probab., 15:no. 18, 526-603, 
2010 . 


30 


[22] L. Erdos, B. Schlein, and H.-T. Yau. Universality of random matrices and local relaxation 
flow. Invent. Math., 185(1):75-119, 2011. 

[23] L. Erdos and H.-T. Yau. Universality of local spectral statistics of random matrices. Bull. 
Amer. Math. Soc. (N.S.), 49(3):377-414, 2012. 

[24] L. Erdos and H.-T. Yau. Gap universality of generalized Wigner and /3-ensembles. J. Eur. 
Math. Soc. (JEMS), 17(8):1927-2036, 2015. 

[25] L. Erdos, H.-T. Yau, and J. Yin. Bulk universality for generalized Wigner matrices. Probab. 
Theory Related Eields, 154(l-2):341-407, 2012. 

[26] J. Huang and B. Landon. Spectral statistics of sparse Erdos-Renyi graph Laplacians, 2015. 
Preprint, arXiv: 1510.06390. 

[27] J. Huang, B. Landon, and H.-T. Yau. Bulk universality of sparse random matrices. J. 
Math. Phys., 56(12):123301, 19, 2015. 

[28] D. Jakobson, S.D. Miller, I. Rivin, and Z. Rudnick. Eigenvalue spacings for regular graphs. 
In Emerging applications of number theory (Minneapolis, MN, 1996), volume 109 of IMA 
Vol. Math. AppL, pages 317-327. Springer, New York, 1999. 

[29] K. Johansson. Universality of the local spacing distribution in certain ensembles of Hermi- 
tian Wigner matrices. Comm. Math. Phys., 215(3):683-705, 2001. 

[30] C.H. Joyner and U. Smilansky. Dyson’s Brownian-motion model for random matrix theory 
- revisited. With an Appendix by Don Zagier, 2015. Preprint, arXiv:1503.06417. 

[31] C.H. Joyner and U. Smilansky. Spectral statistics of Bernoulli matrix ensembles - a random 
walk approach (i), 2015. Preprint, arXiv:1501.04907, to appear in J. Phys. A. 

[32] N.M. Katz and P. Sarnak. Random matrices, Probenius eigenvalues, and monodromy, vol¬ 
ume 45 of American Mathematical Society Colloquium Publications. American Mathemat¬ 
ical Society, Providence, RI, 1999. 

[33] B. Landon and H.-T. Yau. Convergence of local statistics of Dyson Brownian motion, 2015. 
Preprint, arXiv: 1504.03605. 

[34] B.D. McKay. Asymptotics for symmetric 0-1 matrices with prescribed row sums. Ars 
Combin., 19(A): 15-25, 1985. 

[35] M.L. Mehta. Random matrices, volume 142 of Pure and Applied Mathematics (Amsterdam). 
Elsevier/Academic Press, Amsterdam, third edition, 2004. 

[36] S.J. Miller and T. Novikoff. The distribution of the largest nontrivial eigenvalues in families 
of random regular graphs. Experiment. Math., 17(2):231-244, 2008. 

[37] 1. Oren and U. Smilansky. Trace formulas and spectral statistics for discrete Laplacians on 
regular graphs (H). J. Phys. A, 43(22):225205, 13, 2010. 

[38] L. Pastur and M. Shcherbina. Universality of the local eigenvalue statistics for a class of 
unitary invariant random matrix ensembles. J. Statist. Phys., 86(1-2):109-147, 1997. 


31 


[39] M. Reed and B. Simon. Methods of modern mathematical physics. IV. Analysis of operators. 
Academic Press [Harconrt Brace Jovanovich Publishers], New York, 1978. 

[40] M. Shcherbina. Change of variables as a method to study general /3-models: Bulk univer¬ 
sality. J. Math. Phys., 55(4):043504, 23, 2014. 

[41] T. Shcherbina. Universality of the local regime for the block band matrices with a hnite 
number of blocks. J. Stat. Phys., 155(3):466-499, 2014. 

[42] U. Smilansky. Discrete graphs—a paradigm model for quantum chaos. In Chaos, volume 66 
of Prog. Math. Phys., pages 97-124. Birkhauser/Springer, Basel, 2013. 

[43] T. Tao. The asymptotic distribution of a single eigenvalue gap of a Wigner matrix. Probab. 
Theory Related Fields, 157(1-2):81-106, 2013. 

[44] T. Tao and V. Vu. Random matrices: universality of local eigenvalue statistics. Acta Math., 
206(l):127-204, 2011. 

[45] L.V. Tran, V.H. Vu, and K. Wang. Sparse random graphs: eigenvalues and eigenvectors. 
Random Structures Algorithms, 42(1):110-134, 2013. 

[46] B. Valko and B. Virag. Random Schrodinger operators on long boxes, noise explosion and 
the GOE. Trans. Amer. Math. Soc., 366(7):3709-3728, 2014. 

[47] N.C. Wormald. Models of random regular graphs. In Surveys in combinatorics, 1999 (Can¬ 
terbury), volume 267 of London Math. Soc. Lecture Note Ser., pages 239-298. Cambridge 
Univ. Press, Cambridge, 1999. 


32 


